Perturbational approach to the possible quantum capacity of additive 

Gaussian quantum channel 



Xiao-yu Chen 

College of Information and Electronic Engineering, Zhejiang Gongshang University, Hangzhou, 310018, China 



Abstract 



additive Gaussian quantum channel, we have 



For a quantum channel with additive Gaussian quantum 
noise, at the large input energy side, we prove that the one 
shot capacity is achieved by the thermal noise state for all 
Gaussian state inputs. For a general case of n copies in- 
put, we show that up to first order perturbation, any non- 
Gaussian perturbation to the product of identical thermal 
states input has a less quantum information transmission 
rate when the input energy tends to infinitive. 
PACS number(s): 03.67.-a, 42.50.Dv, 89.70,+c 

Quantum capacity is one of the main issues in quan- 
tum information theory. It is concerned with the transmis- 
sion ability of unknown quantum state on a given quan- 
tum channel. The critical quantity involved in the quan- 
tum capacity is the coherent information (CI) I c (a,£) = 
S(£(a)) - S(aQ R ') [I] 0. Here S(g) = -Trplog 2 g is the 
von Neumann entropy, a is the input state, the applica- 
tion of the channel £ results the output state £{<j) ; <j® r = 
(£<S>T)(\tp) with R referred to the 'reference' system[T] 
(the system under process is Q system with annihilation 
and creation operators a and , we denote a® as a for sim- 
plicity), is the purification of a. The quantum channel 
capacity is [3] [4] [5] 



Q = lim sup — I c (a n , £ 

n— >oo _ n 



(1) 



Quantum capacity exhibits a kind of nonadditivity [5] 
that makes it extremely hard to deal with. The first 
example with calculable quantum capacity is quantum 
erasure channel [7]. Other examples are dephasing qubit 
channel [8], amplitude damping qubit channel Sj, and con- 
tinuous variable lossy channel [TU], where the channels are 
either degradable or anti-degradable[ll]. Gaussian quan- 
tum channel [12] (additive classical Gaussian channel fol- 
lowed Holevo [13] ) is quite essential in quantum informa- 
tion theory. Unfortunately, this channel is neither degrad- 
able nor anti-degradable[Hj [E] makes the technics devel- 
oped for calculating the quantum capacity unapplicable. 

The quantum capacity of the Gaussian quantum channel 
has been conjectured as [IB] 



Q = max{0, - log 2 (eiV„)}, 



(2) 



where N n specifies the Gaussian quantum channel. It can 
be achieved by quantum error-correction codes [12]. For 



£&) = ^- I — exp(-\a\ 2 /N n )V(a)crV^a) 

N n J 7T 



(3) 



where T> (a) — exp(aa^ — a*a) is the displacement op- 
erator. Any quantum state a can be equivalently spec- 
ified by its characteristic function Xa-(fJ-) = Tr\aT> (fx)], 
and inversely a = f\Yl i ^ i \Xa{lj)'D{—^)- The charac- 
teristic function of the noisy state a' = £(cr) is Xct(m) = 

A single mode thermal state p has a characteristic func- 



tion of the form x(m) 

p = J^x(vM-p) 



■ ata with v = N/(N + 
v x = N X /(N X + 1)], 



exp[— (N + |) \p\ ], and we have 
(1 - v)v" 

1), [conventionally in the following, 
where N is the average photon number. The noisy state 
is p' = £{p) = (1 — v')v' a a , with average photon num- 
ber N' — N + N n , this accounts for the 'additive' of the 
channel. 

We now consider a Gaussian state pa (which comprises 
thermal noise state as its special case) input to the channel. 
A single mode Gaussian state is described by its real cor- 
relation matrix a (we drop the first moments of the state 
for they can be removed by local operations). We have 



'</</ 
X 'IP 



<-qp 
l PP 



The energy of the Gaussian state is 



E = Tr[((V a + h) pc] = \{o-qq + at pp ). For a Gaussian state 
input pq, the output p' G and the joint output state Pq R 
are still Gaussian. The symplectic eigenvalues [13] of these 
states can be obtained. The coherent information is 

I C ( PG ,£) = g(d - i) - g{di -\)- g(d 2 - \), (4) 



with 



do 
di,2 



^jN n + 2N n E + E 2 x 



(5) 



[N 2 + 2N n E+-±N n D G ], (6) 



where D G = y/{N n + 2E) 2 + 1 - (2E) 2 x, with x = 
det(a)/E 2 . Here g(s) = (s + l)log(.s + 1) - slogs is 
the bosonic entropy function. When E is fixed, the 
maximum of x can be obtained with the derivatives on 
x — \{oL qq + a pp — 2E), where A is the Lagrange multiplier. 
It follows that the maximum value x — 1 is achieved when 



1 



VP 



= y/E, a qp — 0. At sufficiently large input en- 
ergy E, calculating dJc ^' £ ^ and expanding the expression 
with E , we then obtain 



dI c (p G ,£) 
dx 



1 



1 



2Ex 2 "3N n 



2N n ) + o( 



E 2 ' 



(7) 



dx 

is positive for sufficiently large input energy. While x has 
its global maximum value x = 1, so the coherent infor- 
mation achieves its maximum at x = 1 which corresponds 
to thermal noise state input. Hence we can conclude that 
for sufficient large but definite input energy, the one-shot 
quantum information capacity of Gaussian quantum chan- 
nel is achieved by thermal noise state input of all Gaussian 
state inputs. 

For thermal state input p, denote the annihilation and 
creation operators of the 'reference' R system as b and ¥ , 
we have [T(3 

P QR ' = (l-u)(l-v„)exp[^(l-u„)a t 6 t ] 

xv^ a (vv n f b exp[V^(l ~ v n )ab], (8) 

The state p QR can be written as S2(r)(pA <£> Pb)S 2 {t), 
Si(r) = exp[r(a jt b' f — ab)] is the two-mode squeezing 
operator, the squeezing parameter r is determined by 
tanh2r = 2,yN7WTTj/(2N + N n + 1). p A and p B are 
two thermal states with average photon numbers Na and 
Nb, respectively, where Na,b = \{D ± N n — 1) with 
D = ^N 2 + 2{2N + l)N n + 1. The coherent information 
will be [13! [IS] 



I c (p,£) = g(N + N n ) - g(N A ) - g(N B ). 



(9) 



One of the useful formula is 

N = 7V B cosh 2 r+ (N A + l)sinh 2 r, (10) 

To treat with the non-Gaussian perturbation, we need 
the following lemmas: 

Lemma 1: (£ l){a) k p QR a m )=v -(k+m)/2 b k p QR' ¥ m^ 
proof: With the characteristic function \ oi p^ R , the 
lhs can be written as / exp[-J^- - 

(at - a*) k x QR {p) V(~p)(a ~ a) m . After the integral on 
fi = (/ii,/X2) 5 and a displacement on a (in the ordered 
operator product): a — ► a + a; a* — ► a* + a* , the lhs can 
be further written as 



1 

An 



d 2 c 



tfc exp[-^ + a( — + 



• exp^-^ + v^-^-^K 



where the notation : H : refers to that H is in its ordered 
operator product, that is, all creation operators are at the 
left of the annihilation operators. The integral can be 



worked out with the formula / = -k J exp[— +aa + 
a*r] = exp(Kar), and its derivatives 

flfc+mr mi„{^ m } 



da m dr k ~ ^ \U\l 

1=0 



Since N n < 1/e (Otherwise I c = 0, this is shown by We have 

also by numeric results of Gaussian input), thus dJc( - PG ' £ - ) 



(12) 



min{/c,m} 

(£®I){<J k p QR a m )=: ^ 



1=0 



Or + ^b) k -\ir + y/ttfr-'pQ* : . (13) 



1 A n 



Removing the ordered notation by moving all creation op- 
erators to the left of p® R and all annihilation operators 
to the right of p® R , and exchanging b with pQ R according 
to 

(fl*b = [b- ViJ(l - v n )at]/(vv n )(fl* , (14) 

exchanging with 6 and further with p® R , where the 

formula tf m b n = E^ {m ' n} (?) (?)*!(-l) < 6 n_< & t(m ~ <) is 
used, after all the summation, the lemma 1 is proved. 

Lemma 2: {£ ®l)(a k p QR a^ m )^v l ~ k+m ^ 2 b^ k p QR ' b m . 

proof: The lhs is £ / ^(a - a) fe {/ ^ exp[-^ + 
Ml a*-a^] x ^( / x)©(- Al )}( a t- ( 
on /z, it is 



l , after the integral 



d 2 a 



■exp 



-](a-a) fc 



1 

A^,/ 7T "~" L A„ 
expftW + y^a* - a*)&t] |00) (00| 
exp[a*a + ^/v(a - a)6](a t - a*) m , 



(15) 



where the formula : e a a := |0) (0| is used. Note that 
(a — a) e aa = e aa a, thus (a — a) fc exp[aa t + ^/v(a^ — 
a*p] joO) = exp[aat + ^(ot - a*)&t](^fet)fc | 0) . After 
the integral on a, the lemma 2 is proved. 

In the single mode situation, we expand the input state 
p s at the vicinity of p, the characteristic function of the in- 
put state is xM =Tr(p eJ D(^)) = x(p)Q-+ef(j*, M*))- The 
perturbation item f(p,p*) is a polynomial of /x and /x*. 
Typically, this may contain (1) \p,\ 2n (n > 1, the n = 1 is a 
Gaussian type perturbation, here we discuss single mode 
situation, thus n can not be confused with that appeared 
in (fT]) where it stands for the number of the modes ) and 
(2) cp n {-p*) l +c*p l (-p)* n (n ^ I). The first type of per- 
turbation will contribute to the first order perturbation of 
the eigenvalues of p, while the second type of perturbation 
has not a first order perturbation to the eigenvalues of p, 
it can only contribute to the second order perturbation. 
For the first type perturbation \p,\ 2n , We have Xe{p) = 
(11) x(M)(l + s\p\ 2n ), thus p E = Jexp[-(N + |) \p\ 2 }(l + 
e\p\ 2n )V{-p)^ = (l +£ (-ir^)p = + 8^. The 



strict eigenvalues of p e are A 



k 



)p = p + eq 

Afc + £0fe, with 
Afc = (1 - v)v k and qb k = A fc ^ fc ,where £ fc = (1 - 
v) n T,^ {n ' k} (-l) j nl^) Q)N-J. £t is the eigenvalue of an 



2 



operator 



£(a,a + ) = (l-i))*3-l) 3 *^ "].V ' f / T; ./ ; - 



(16) 



with its eigenvector |fc) . The moments of <p can be cal- 
culated by Tr{a) l a l cj>) = (-l) n ^Tr(a^ l a l p) = (-1)" 
-Ar^(l\N l ), thus for / < n, the moments are nullified. The 
entropy of p e can be expanded up to the second order of e 
as S(p s ) = S{ P ) - \e 2 E fe ^/A fe + o(e 3 ), where Trfa) = 
and Tr(a<a4>) — are used. We have Ylk&k/^k = 
Trm = (l-«) w E?=o(-l)^(i)^~ J Tr(otV0) - (1- 

v) n N~ n (n!) 2 . The calculation of the entropy of the noisy 
state is straightforward, it is 



S(p's) = S(p') 



(n\y 



1 



2 iV" l (7V' + l)' 



+ ( £ 3 ). 



(17) 



The purification of p is p® R — J2km V AfcA m \kk) (mm\ . 
The state pf R then is expanded in e to the linear item 
(the e 2 term is less important in the large iV limit). 
p QR = pQR + with $ = i( $o + $ t) 5 $ Q = (! _ 

«)™E" = o(- 1 ) J 7f(™) Ar ^ ab&tJ ^ /2 P Qfl - Using lemma 1, 
we arrive at 



(£ ® I)p? fl = p QR ' + e$', 



(18) 



with $' = |($' + % = £,(b,tf)p QR ' . The eigen- 

states of p® R are |fcm) = 5*2 (r) \ km) with eigenvalues 
A/cm = (1 — u a)i^ (1 — vb)vb - The corresponding anni- 
hilation operators are transformed to a' — S^r^aS^r) = 
a cosh r — 6^sinhr, b' = S2{r)bS 2 {r) = b cosh r — a^sinhr. 
The first order perturbation to the eigenvalues will be 
Km = (km\'&\km)' = \ km {kmU{b,tf)\km)'. Note 
that 



(km\' b^V \km)' = (fem|'£ V- ) b'H'^'^-^a'^^ 



cosh 21 r sinh 2 ^ 8 ' r Ifcr 



(19) 



We can construct an operator SI which is the diagonal part 



of in the basis \kr 



Q = N- n (N + iy 



that is, 



[A^(A^ + l)sinh 2 r\ 



[N B (N B 

d" 



=7P 



dN B dNl~ J 



1) cosh 2 rp 
(20) 



Where we have used the following lemma. 

Lemma 3: 1. a^Vp = N^L)^i N + ^IW' 2 - 

aiaVp=(N + lYYLo%(?)X i & 

Proof: These two equalities can be proved with mathe- 
matical induction. 

For pQ R is a direct product state in the basis \km) , we 
can treat the two new modes separately. Up to e 2 item, the 
entropy will be S(pf#) = S(pQ R ')-ie 2 (£ fcm <&L/A fcm + 



77) + o(e 3 ), where rj is the item comes from the e 2 item in 
the expansion of pf R on p® R , this item can be omitted 
comparing with the main item J2km ^'km/^km f° r large N. 
Denote flkm — (km\ SI \km) , and make use of Eq.(f5D|), we 
have 



5>' fc 2 m /A fcm = $>L/A 



km 



km 



km 



(n\) 2 N- 2n (N + l)- 2n 



E 

j=0 



2j 



(21) 



where Eq. ([TU|) has been used, and B — Nb(Nb + 1) cosh 4 r, 
A = N A (N A + 1) sinh 4 r. At the limit of N — » oo,we have 
iV' — * iV, the difference of CI between p £ and p is 



lim [/ c (p 6 ,£)-/ c (p,f)] 

iv— »oo 



-e 2 (nl) 2 N- 2n 
2 v ; 



n 

x[l-2~ 2n Y, 

3=0 



< 0. 



(22) 



If we have a linear combination of above type perturba- 
tions, see ci \p\ 2ni + C2 \p\ 2n2 (n-i < 112), the interference 
term J2 k fakfah I % = (-l)™ 2 ^^^) - 0, the similar 
zero interference can be proved for the perturbation to the 
joint state p® R . Thus each item act separately to the co- 
herent information. So that at the sufficiently large input 
energy and at the single use of the channel, the first or- 
der non-Gaussian perturbation to the input thermal state 
will not improve the conjectured capacity of the Gaussian 
quantum channel. 

In the n use of the channel with an input Gaussian 
state p n ,the algebraic equations of the symplectic eigenval- 
ues [13] are not analytically solvable. Fortunately, when 
the input state is a two mode squeezed thermal state, we 
can obtain analytical result. For the input state of p 2 = 
£2(^2) P® pS\{r2) , where S , 2(^2) is the two mode squeezing 
operator with real parameter T2 , the coherent information 
is I c (p 2 , £ ® 2 ) = 2 max{0, g(d - 1) -g{di - %)-g(<k - 5)}, 
where in the expression of dj, E should be substituted with 
E/2, now E is the total energy of the two mode state; and 
x = 1/ cosh 2 ?'2 . The detail of the calculation will be given 
elsewhere. The conclusion is that the maximum coherent 
information is achieved by the product thermal state p®p 
when the input state is two mode squeezed thermal state 
at sufficiently large input energy. Also with Lagrange mul- 
tiplier method, it is easy to prove that at sufficiently large 
input energy, for all n mode product thermal states, the 
product of identical thermal state will achieve the maxi- 
mal of coherent information. Thus an unbalanced energy 
distribution among modes will not increase the total co- 
herent information. Hence, we have proved that for all 
product Gaussian state inputs and all product of two mode 
squeezed thermal state inputs, the maximum of the coher- 
ent information is achieved by product identical thermal 
state p n — p® n for sufficient large input energy. 

We now turn to the first order multi-mode perturbation. 
We will omit the case that the perturbation can be treated 



3 



separately for each mode. What left is the perturbation 
with particle exchange among modes. A typical case is 

Xne(p) = + s(c^ ^ ■ ■ ■ flfrtf 1 ftf ■ • ■ Mn" + 

C*l4 kl l4 k * ■ ■ ■ /»n fen /»lM2 ' ' ' )] with EiLl 

zC"=i k — m (the requirement of first order pertur- 
bation). We may denote the perturbation as (k, 1), with 
vectors k = (ki, k 2 , ■ ■ ■ , k n ), 1 = (h,h, ■ • ■ , In) and k^l . 
The perturbed input state is p ne = p® n + (f>. To simplify 
the calculation, we introduce a generation function 



U (r, a) 
then 



n, 



X„(ju)D(-pt)exp[/i- 



d 2 "V (r, g) 

n^d^) 



9 2m / (r, a) 



' / -! /T— fl 



(23) 



The eigenspace of p ne can be classified as subspaces ac- 
cording to the eigenvalue p® n . When the eigenvalue of 
p®" is Xj = (1 — v) n yi, the basis of the subspace can be 
\j1J2, ■ ■ ■ ,j n ) = |j) j with Yh=iM = 3- Tlie action of <f> 
will keep j invariant, that is, it is an operator in j subspace. 
In the subspace, we may specify 4> as Mj. The eigenvalues 
of Mj is supposed to be Aji, then £\ A ji = Tr(Mf) = 

EjOlf|j)ai0|j> =E J j'<il0|j')a , ^|j> = Ej(j|0 2 |j>, 

where j' may not be in the j subspace. Here we have used 
the fact that (j| <f> = for j' 0he subspace of j. The 
perturbation to the entropy will be 

S(p ne )-S(p®") = -l £ iJ2^+o(e 3 ) 

= -\e 2 Tr(cf > 2 /p^)+o(e 3 ).(24) 

Where the linear term of e is nullified by the fact that 
Tr(Mj) = 0. For 



Tr[I <t> (T,a)It(T , ,a')/p® n ]=e X j>[- 



n N + r 



(25) 



Thus 



Tr((J) 2 /p® n ) = 2 |c| 



2 iWZj!) 
[N{N + 1)Y 



In obtain Eq. ([25| . we first work out the operator integral 
of / part, which have an operator that conceal the \j p® n 
operator. Eq. ([25|) exhibits that any interference item of 
Tr(W j p® n ) type will be nullified for (k, 1) ^ (k', 1'). Thus 
each perturbation item contributes to the entropy sepa- 
rately. 

The perturbation to the joint QR state is more sophisti- 
cated. With almost the same routine as we do in obtaining 
the perturbation operator in one mode situation, we can 
get p® R = pQ R ® n + £$. Care should be taken in obtaining 
the operator $ for complex coefficient c, when c is real, 
there are no problem; while c is complex, there are ex- 
tra phase factors in transforming the eigenbasis of <f> into 
the basis of direct product of the modes. However the 
phase factor can be absorbed in the purification process, 
that is, if we have e l8 |j) |j) in the purified state, it will 
make no difference with |j) |j) in obtaining the reduced 



state. We thus have $ = -|($o + $0), The generation 



2 

exp 



function of $0 is I<s> 
The action of the channel then is 1'^ 



N+l 



exp [^S] 

= {£ ® I)U = 
exp(pr • b^) exp[-^j- — pa ■ h]pQ R ®" according to lemma 

1 and lemma2, where p = [N(N + l)] -1 / 2 . The contribu- 
tion to the entropy should be evaluated in the eigenbasis of 
pQR yy e ma y denote the subspace of p QR ® n as|i, i; j, j)' 
which has eigenvalue = (1 — i;^)"(l — VB) n v\v B . In 
this subspace, we denote <&q as My, the elements of 
are (i, j| $q \i, V; j, j') , the sum of the square of the 
eigenvalue of is TrM 2 -. We obtain the contribution to 
the entropy by first evaluating (i, j|' 1^ \i, j')' = 

(i, i-j, j| v^ n r^v® n I*. >';iJ') > which is ex p[ T ■ °/( N + 1 )] 

(i, j| exp[pr • (b^ coshr + asinhr)] exp[— pa ■ (bcoshr + 
a* sinhr) ■ We expand the exponent of the op- 

erators to drop the terms that do not keep the total 
particle numbers in A and B parts respectively. De- 

note I B {r,a) = £~ ^(pcosh r) 2fe (r • h f ) k (a ■ b) fc , 

Ia(t, a) = Efco ^(psinhr) 2fc (r • a.) k (a -^) k . Then in 
the calculation of the contribution to the entropy become 
a trace on the whole space, the restriction on the subspace 
is removed. We have the generation function 



F 



exp[(T • a + t / -o / )/(JV + 1)] 
xTr[(I B (r,a)I B (r',a')/pT] 
xTr[I A (r,a)I A (r',a')/pT} 
00 m B l A m ~ l 

m=0 i=0 



[l\{m-iy.] 2 [N(N + 



(27) 



, with which the perturbation to the entropy can be calcu- 
lated, where we have used another generation function in 
evaluating F, and at the final step we exchange the orders 
of summation and make use of Eq. (fT0|) to conceal the fac- 
tor exp[(r • a + t' -a') / (N + 1)]. The generation function is 
about that the two ingredients both come from $0 , if both 
(26) come from $J,the result will be the same. In the case inter- 



cross of $0 and $J, we should substitute (r • a') m (r' -a) 1 
inEq.lgll) by {t ■ T r ) m {a' -a) r 



Tr[\{M l3 + 7l4)] 2 



= 2 Id 



[A^(A^ + 1)] 2 



X J2[ ) S'^ m_ '(28) 

The situation is strictly like that of the single mode case. 
Eq. (|2T|) indicates that each perturbation term contributes 
to the entropy separately. Also, the state pQ R ® n then 
is expanded in e to the linear item (the e 2 term is less 
important in the large N limit). 

We have shown that all first order perturbation to the 
input product identical thermal state can only decrease the 
coherent information at large input energy. In one mode 



4 



case of Gaussian input, in two mode of squeezed thermal 
input, and non-Gaussian perturbation of one mode as well 
multi-mode to thermal state input, thermal or their identi- 
cal product achieve maximum of the coherent information 
at infinitive input energy. For all inputs of product Gaus- 
sian states, product two-mode squeezed thermal states, we 
have proved that the channel capacity of additive Gaussian 
quantum channel is described by conjectured formula (J2|) 
and achieved by product of identical thermal state. All 
kinds of the first order perturbations (non-Gaussian) to 
the product of identical thermal state input will lead to a 
less coherent information for sufficient large input energy. 
Whether a bigger coherent information can be achieved by 
a state which is far from the product of identical thermal 
state still remains open. 
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